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In this article we compute the density of Dirac particles created by a cosmological 
anisotropic Bianchi I universe in the presence of a constant electric field. We show that 
the particle distribution becomes thermal when one neglects the electric interaction. 

1. Introduction 

During the last decades a great effort has been made in understanding quantum 
processes in strong fields. Quantum field theory in the presence of strong fields 
is in general a theory associated with unstable vacua. The vacuum instability 
leads to many interesting features, among them particle creation is perhaps the 
most interesting non-perturbative phenomenon. After the publication of the articles 
by Parker^and Zeldovich^, the study of quantum effects in cosmology became 
a very active research field. In the last thirty years a large body of papers has 
been published on the problem of vacuum effects in isotropic and homogeneous 
gravitational backgrounds, mainly in deSitter and Robertson Walker models, and 
only a few try to discuss quantum processes in anisotropic Universes. 

In order to analyze the mechanism of particle creation in cosmological back¬ 
grounds we have at our disposal different techniques such as the adiabatic approach^ 
the Feynman path integral method ^, the Hamiltonian diagonalization technique®d^ 
as well as the Green function approach® 

The study of quantum effects in gravitational backgrounds with initial singular¬ 
ities presents an additional difficulty. The techniques commonly applied in order 
to define particle states are based on the existence of a timelike Killing vector or 
an asymptotically static metric®. A different approach is needed to circumvent the 
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problem related to the initial singularity. In this direction, the Feynman path- 
integral method has been applied to the quantization of a scalar field moving in 
the Chitre-Hartle Universe®’^^. This model has a curvature singularity at t = 0, 
and it is perhaps the best known example where a time singularity appears and 
consequently any adiabatic prescription in order to define particle states fails. A 
spin 1/2 extension has been considered by Sahni^^. 

A different approach to the problem of classifying single particle states on curved 
spaces, is based on the idea of diagonalizing the Hamiltonian. This technique per¬ 
mits one to compute the mean number of particles produced by a singular cosmo¬ 
logical model, and in particular by the Chitre-Hartle Universe®. 

An interesting scenario for discussing particle creation processes is the anisotropic 
universe associated with the metric 

ds^ = —dt^ + t^^dx^ + dy"^) + dz^ ( 1 ) 

The line element (1) presents a space-like singularity at t = 0. The scalar curva¬ 
ture is R = 2/t^, and consequently, the adiabatic approach® cannot be applied 
in order to define particle states. Creation of scalar particles in the anisotropic 
universe (1) was originally computed by Nariai with the help of Feynman propaga¬ 
tors, obtaining as a result that the creation process occurs in accordance with the 
black-body thermal law in a 2-dimensional hypersurface related to the anisotropic 
cosmic expansion. With the help of the Hamiltonian diagonalization method®’^®’^'^, 
Bukhbinder^ could compute the rate of scalar particles produced in the space with 
the metric (1), obtaining as a result, a Bose-Einstein distribution. In Ref.^® this re¬ 
sult has been extended including a time dependent electric field More recently^®’^^, 
a quasiclassical approach has been applied to compute the rate of scalar as well as 
Dirac particles in the cosmological universe associated with the metric (I). 

The introduction of an external electric field permits one to consider an addi¬ 
tional source of quantum processes. The density of particles created by an intense 
electric field was first calculated by Schwinger^®, different authors®’^® have discussed 
this problem. Pair creation of scalar particles by a constant electric field in a 2-1-1 de 
Sitter cosmological universe has been analyzed by Garriga®®. Quantum effects asso¬ 
ciated with scalar and spinor particles in a quasi-Euclidean cosmological model with 
a constant electric field are discussed by Bukhbinder and Odintsov®^. It is worth 
mentioning that the presence of primordial electric fields could enhance the particle 
creation mechanism and also produces deviations from the thermal spectrum. 

The purpose of the present article is to discuss the production of Dirac particles 
in the anisotropic cosmological background associated with the line element (1) in 
the presence of a constant electric field. In order to compute the rate of particles 
created we apply a quasiclassical approach that has been successfully applied in 
different scenarios^"^’®®’®®’®^. The idea behind the method is the following: First, 
we solve the relativistic Hamilton-Jacobi equation and, looking at its solutions, 
we identify positive and negative frequency modes. Second, we solve the Dirac 
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equation and, after comparing with the results obtained for the quasiclassical limit, 
we identify the positive and negative frequency states. 

The paper is structured as follows. In Sec. 2 we solve the relativistic Hamilton- 
Jacobi equation and compute the quasiclassical energy modes. In Sec. 3, after 
separating variables, we solve the Dirac equation and obtain the density of Dirac 
particles created. The discussion of the results and final remarks are presented in 
Sec. 4. 

2. Hamilton-Jacobi equation 

The relativistic Hamilton-Jacobi equation can be written as 

where here and elsewhere we adopt the convention c = 1 and h = \. 

The vector potential Aa 


Aa = (0,0,0,-Et), (3) 

corresponds to a constant electric field Ek. The corresponding invariants = 

—2E‘^ and F^'^*F^^ = 0 indicate that there is no magnetic field. Since the metric 
gaf) associated with the line element (1) only depends on f, the function S can be 
separated as 

S = F{t) + kxX + kyy + kzZ. (4) 

Substituting (4) into (2) we obtain 

^ + {kz + eEt)"^ + rri^. (5) 

The solution of Eq. (5) presents the following asymptotic behavior: 

lim F = ±—t\/e'^E'^t'^ — ^ log{eEt + \/ e^E'^E — w?), (6) 

i^oo 2 2eE 

$ = ^ (eEt)^^ (7) 

as t ^ oo, and 

liin F = ±^(fc2 + fc2)logt, $ = ^ (g) 

as t ^ 0, that is, in the initial singularity. Notice that the time dependence of the 
relativistic wave function is obtained via the exponential operation $ ^ exp(fS'). 
Here the behavior of positive and negative frequency states is selected depending on 
the sign of the operator idt- Positive frequency modes will have positive eigenvalues 
and for negative frequency states we will have negative eigenvalues. Then in Eqs. 
(6), (7) and (8), upper signs are associated with negative frequency values and the 
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lower signs correspond to positive frequency states. After making this identifica¬ 
tion we proceed to analyze the solutions of the Dirac equation in the cosmological 
background ( 1 ). 

3. Solution of the Dirac Equation 

The covariant Dirac equation in curved space in the presence of electromagnetic 
fields can be written as follows 

+ = 0, (9) 

where the curved gamma 7 ^ matrices satisfy the anticommutation relation { 7 ^^, 7 ^} = 
2g^'' and the spinor connections are 

= (i») 

where 

s"^ = ^(7V-7V). (11) 

The matrices establish the connection between the curved 7 ^ and Minkowski 

7 ^ Dirac matrices as follows 


7m = 7^ = aj)7^ 

(12) 

7 ^ 7 " + 7"7^ = 

(13) 


Since the line element (1) is diagonal, we choose to work in the diagonal tetrad 

7 ^ = \/|ff^''l 7 ^) no sum. (14) 

Substituting (14) into (10), we obtain 

ri = i7V, r2 = iff, r3 = o, ro = o (i5) 

and the Dirac equation (9) takes the form 

{’“I +s + +“}«'» = ». (“) 

where 'I'o = The factor t was introduced in order to cancel the contribution due 
to the spinor connections (15). The equation (16) can be written as a sum of two 
first order commuting differential operators as follows^®’^® 

(Xi -f k2)<^ = 0 

= M> = -ki<^>, 


(17) 

(18) 
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where the spinor $ is related to 'I'o via the equation 


and A: is a separation constant. The operators Ki and K 2 read 

d d 

Ki = t 7°(^ - ieEt) + m7^7° 


(19) 


( 20 ) 


Since Eq. (16) commutes with —iV, the spinor $ can be written as $ = <I)o exp(z(A:xa:+ 
kyU + kzz)). Choosing to work in the following representation of the Dirac matrices 



( 22 ) 


we obtain that K 2 ^ = reduces to an algebraic equation that permits one deter¬ 
mine the relation between the components of the bispinor $ 


$0 — 


$1 

$2 



$1 

kx(7^ ( 
iky — k 


(23) 


where eigenvalue k is 

k=ifkl + (24) 

Using the representation (22) and taking into account the spinor structure (23) we 
obtain that, for kz = 0, Eq. (20) reduces to the system of equations 


dipi 

dt 


+ -t/’i + {m + ieEt)(p2 = 0, 


(25) 


d(fi2 k , , , , 

—— -I- -ip 2 + (m — ieEt)ipi = 0. (26) 

In this way, we have reduced the problem of solving Eq. (16) to that of finding 
the solution of Eqs. (25)-(26). From (25) and (26) we obtain the following second 
order differential equation 


f]h.f 

dt^ 




2fc-b3/4 

I 2 



V'2 = 0, 


(27) 


where we have introduced the variable ■i /'2 = and have neglected the mass 

in the first-order variation of (j) 2 - The solution of (27) can be expressed in terms of 
Whittaker functions as 


V'2 = CiMx,y{ieEf) + C2Wx^f_t{ieEf), 


( 28 ) 
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where 


im? k 1 

“~ 4 ^’ 2 ~ 2 


( 29 ) 


and Cl and C 2 are arbitrary constants. Analogously, after introducing the new 
variable ipi = we have that the equation for ipi is 


dt^ 


k^ + 2k + 3/4: 


) t/-! = 0 


(30) 


with solutions 

V’l = C3 Ma.^+i (teAt^) + CiWx,^+i{ieEf) (31) 

where the Whittaker functions Wk,ij.{z) and Mk,fi{z) can be expressed in terms of 
confluent hypergeometric fuctions^^ as follows 


MkA^) = +^i-k,l + 2^J^, z), (32) 

WkA^) = e-^/2zi/2+M|y( 1 + ^ 1 + 2^^, z). (33) 

In order to construct the positive and negative frequency modes we use the asymp¬ 
totic behavior of the hypergeometric functions (32) and (33) and compare the so¬ 
lutions of Eq. (27) with (6) and (8), obtained after solving the Hamilton-Jacobi 
equation. Using this procedure and looking at the asymptotic behavior of Mk,fj,{z) 
as z ^ 0 

Mfe,^(z)-e-^/"zi/2+M (34) 

and using the fact that all the coefficients in Eq. (27) are real, we obtain that the 
positive and negative frequency solutions as t ^ 0 are 


(35) 

where is a normalization constant. Analogously, looking at the behavior of 


Wk,p,{z) as |z| ^ 00 

lAfe,4z) - e-^/V, (36) 

we have that the corresponding positive and negative frequency modes as t ^ +oo 
are 

= €toWkA^eEf )^ = £-VE_fc,^(-ze£;t2), (37) 

where and are normalization constants. 

Using the asymptotic expressions 

lim M(a, 6, —z/a)/r(5) = z 2“2 J{,_i(2-\/z) (38) 

a —>-oo 

lim U{a, b,-z/a)T{l + a - b) = z^-^Hfjf (2v^) (39) 

a^oo 
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we obtain that, as E ^ 0, the positive and negative frequency modes (35) and (37) 
reduce to those obtained through using the diagonalization method ^ as well as with 
the help of the semiclassical approach^® 

The positive frequency mode as t ^ 0. can be expressed in terms of the positive 
and negative frequency modes via the Bogoliubov transformation 

V'o" = aV’i + (40) 

The Whittaker function Mfc^^(z) can be expressed in terms of Wk,ii{z) as^"^ 

(41) 

Using the expression (41), we have that the negative frequency solution tp^ can be 
written in terms of and ■0^ as follows 


r(2M + i) , r(2Ai + i) 


r(/r - A + 4) 




r(i + ^ + A) 


( 42 ) 


where we have made use of the property W-k,fi{—z) = (M4^^(z))* 

Since we have been able to obtain single particle states for in the vicinity of 
t = 0 as well as in the asymptote t ^ oo, we can compute the density of particles 
created by the gravitational field. With the help of the Bogoliubov coefficients®’®. 
From (42) and using the fact that tpQ = q; 0“ + /3(0“ )*, we obtain. 


_ g2'i7rAi 


r(4 + ^ — A) 
r(4 + ^ + A) 


Substituting into (43) the values for /i and k we obtain 

l«l" (^ + ii)sinh(^ + f^)’ 


where we have used the relation^® 


|r(*2/)|" = 


y sinh Try 


The computation of the density of particles created is straightforward from (44) 
and the normalization condition®® of the wave function 
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It is worth mentioning that, thanks to the relations (43) and (46) we did not 
have to compute the normalization constants and 


4. Results and Discussion 

Let us analyze the asymptotic behavior of (47) when the electric field vanishes. 
Taking into account that sinh{z) ~ /2 as 2 ^ oo, we readily obtain 


n « = exp(-27rfc_L), (48) 

|a| 

which is the result obtained in^®. Expression (48) corresponds to a two dimensional 
Fermi-Dirac thermal distribution. In the case of strong electric fields the density 
number of scalar particles created takes the form 


n 


^ — Tzk\ — 


(49) 


showing that the density of particles created by the cosmological background and 
the electric field (49) is a Fermi distribution with a chemical potential proportional 
to Integrating the particle density n (49) on momentum we obtain the total 
number of particles created per unit volume. 



ndkxdkydkz 


1 

t2(27r)2 


nk±dk±dkz. 


(50) 


In order to carry out the integration we have to notice that n does not depend 
on kz and consequently integration over kz is equivalent to the substitution®’^® 
f dkz eET, where T is the time of interaction of the external field. Substituting 
(49) into (50) , we obtain that the total number N of particles per unit volume 
takes the form 

Result (51) resembles the number of particles created by a constant electric field 
in a Minkowski space®’^®. It is worth mentioning that the number N of particles per 
unit volume is inversely proportional to T~^ and vanishes as T —> 00 . The volume 
expansion of the anisotropic universe (1) is faster than the particle creation process, 
therefore N becomes negligible for large values of T. Since the gravitational density 
p = l/(87rt2), decreases faster than (51), the particle creation mechanism effectively 
isotropizes in the anisotropic universe (1) in the presence of strong electric fields. 

The results (49), and (51) show that the anisotropic cosmological background 
(I), as well the constant electric field, contribute to the creation of scalar particles. 
The quasiclassical method gives a recipe for obtaining the positive and negative 
frequency modes even when spacetime is not static and an external source is present. 
The presence of the anisotropy with a constant electric field gives rise to a particle 
distribution that is thermal only in the asymptotic field regime. The method and 
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results presented in this paper could be of help to discuss quantum effects in more 
realistic anisotropic cosmological scenarios. 
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